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INVESTIGATION OF THE PERFORMANCE CHARACTERISTICS OF A 
DOPPLER RADAR TECHNIQUE FOR AIRCRAFT COLLISION HAZARD WARNING 


I. SUMMARY 

This report presents the results of a computer simulation study of a 
cooperative doppler radar system for aircraft collision-hazard warning 
under development by NASA-Langley Research Center personnel. The system 
performance under multiple aircraft conditions, with particular emphasis 
on the problems of saturation and interference, is evaluated by simulation 
of the system using a data base consisting of twelve one-hour samples of 
aircraft position data obtained from the ground radar at the Atlanta, 

Georgia airport. In addition, the system error sources are determined by 
comparison of simulated data with flight test data, and other miscellaneous 
system problems are considered. Statistical studies of the data base have 
been previously documented in the Phase I final report [1] and Phase II 
interim report [2] of this contract. 

II . INTRODUCTION 

In phase I of this contract, the feasibility of a dynamic computer 
simulation of the aircraft collision hazard warning system developed by 
NASA-Langley personnel was investigated. The simulation was found to be 
feasible and well within the capabilities of the available digital computers. 

One of the major objectives of the phase I effort was to collect data 
defining aircraft motions in a typical airport terminal area. With FAA 
cooperation, twelve one-hour samples of digitized radar data, controller- 
aircraft voice communications tapes, and controller log sheets were collected 
at the Atlanta, Ga, terminal. These data were taken during morning, 
afternoon, and evening peak traffic periods over a five day interval 
during the month of August 1967. The data were edited by FAA personnel at 
National Aviation Facilities Experimental Center (NAFEC) and made available 
to RTI in the form of digital magnetic tape. The edited data contain, for 



the majority of the aircraft within a 35 mile radius of the Atlanta Airport, 

(1) position data in xyz coordinates for all aircraft under track at four 
second intervals (approximately 700 aircraft tracks with a total flight time 
of approximately 119 hours), (2) coordinate rates at four second intervals 
for all aircraft under track and, (3) supplementary data such as 
identification of aircraft and time of day. The data base is described in 
detail in the Phase II Interim Report on this contract [2]. 

The objectives of the Phase II efforts have been to: (1) development 

of mathematical models of the collision avoidance system and development of 
computer simulation programs to evaluate the performance of the system; 

(2) to use the data describing a typical terminal area to conduct a detailed 
statistical analysis of the effectiveness of various parameters (i.e., 
normalized doppler rate, time to closest approach, or projected miss 
distance) in reducing the probability of false alarms; and, (3) to provide 
assistance in the planning and evaluation of flight test data by use of the 
simulation model. 

The effort described under (2) above has been reported upon separately 
in the Phase II Interim Report. Hence, this report presents the studies 
having the objectives (1) and (3) described above. 

Two experimental collision warning systems have been constructed by 
LRC personnel and have undergone flight test evaluation. This evaluation 
provided a good indication of the system performance when only two aircraft 
are involved. The computer simulations described in this report evaluate 
the system performance under multiple aircraft conditions and provide 
quantitative results on the severity of the problem of saturation due to 
multiple aircraft interrogations. 

In the following, the simulation techniques used are described in 
detail and the analytical models documented. The results of the experimental 
flight tests are compared to simulated results and the sources of system 
error are determined. The severity of the problem of system saturation and 
interference in a multiple aircraft environment is investigated in detail, 
and conclusions and recommendations are given for future work. 
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Ill . SIMULATION TECHNIQUES 


A* SYSTEM DESCRIPTION 

Figure 1 indicates the basic operation of the system. The protected 
aircraft is equipped with a radar unit consisting of two transmitters, a 
receiver , and a display console, and all cooperating aircraft have a 
transponder. The protected aircraft transmits a pair of signals at 
different frequencies. These signals are received at the intruding 
aircraft and multiplied to obtain a difference frequency. The difference 
signal is then re-transmitted to the protected aircraft where its frequency 
is compared with the difference of the transmitted frequencies to obtain 
the Doppler frequency as a measure of closing velocity. Since the 
transponder output power is proportional to the product of the power of 
the two interrogating input signals * the signal power returned to the 
receiver varies inversely with range to the sixth power. This relatively 
sharp fall-off of received power with range permits a relatively accurate 
measure of range to be obtained from measurement of returned signal power. 

To provide for multiple access to the transponder, pseudo-random 
transmitting coding is used. Both of the transmitted signals are randomly 
frequency modulated; a single frequency is derived at the transponder by 
multiplication and filtering, and the random modulated return signal is 
decoded at the receiver. The signals due to multiple aircraft interroga- 
tions appear at the ownship receiver to have a noise- like characteristic, 
while the ownship returned signals are essentially stationary in frequency. 

Because of the range limitation of the radar and the use of random 
transmitter coding, all radars may operate on the same frequencies, 
eliminating the need for separate channel assignments for each aircraft. 
Also, because of the C.W. mode of operation, the peak transmitter powers 
are relatively low, allowing the use of solid-state equipment, and there 
is no specific requirement for a high degree of accuracy in any of the 
frequencies involved . 

This system is primarily a cooperative, C.W. Doppler radar which 
provides a precise measure of the closing velocity of an intruding aircraft, 
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and a somewhat less precise measure of the relative range between aircraft* 
The relative range divided by the closing velocity is computed, and gives 
a measure of the "time to collision." An indication to the pilot of the 
direction to a hazardous target is provided by a multiple lobe receiving 
antenna and associated signal processing to provide nine separate indicators 
of relative target bearing (e.g. up-right, up-center, ahead-right, 
down-left, etc.). 

Figure 2 shows a more detailed block diagram of the system. Not shown 
on this diagram are provisions for filtering the ownship transmitter signal 
from the ownship transponder and the ownship transponder return from the 
ownship receiver. The operation of the system under multiple interrogations 
may be understood by inspecting the spectral densities in Fig. 3. The 
transmitter outputs are composed of the sum and difference of the two 
frequency modulated oscillators. The transponder receives a pair of input 
signals from every radar interrogating it, and will generate an output 
signal for every cross product that exists. For instance, if two radars 
are interrogating the transponder, four cross products will be generated. 
Only two of these products are desirable, and to suppress the undesired 
products, the transmitted signals are randomly frequency modulated such 
that their power is evenly distributed over the transmitter bandwidth. This 
modulation is applied such that it is coherent between signals in a pair 
and not coherent between signals not in a pair. Thus, the noncoherent 
products are spread over a band roughly two times that of the transmitter 
bandwidth, whereas the desired returns from the transponder are spread 
over a band determined by the deviation of FM oscillator no. 2. A filter in 
the transponder passes only the desired signals and the noncoherent power 
contained within the desired signal bandwidth. 

At the receiver, the incoming signals are mixed with the signal from 
FM oscillator no. 2 to derive the desired returns. The desired returns 
are thus stationary in frequency at the Doppler filter output (except for 
the Doppler shift). Signals entering the receiver which are returns to 
other receivers in the area are noncoherent with the mixer signal and are 
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Figure 1. Functional system diagram* 
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Figure 2. Simplified block diagram of transmitter, 
transponder and receiver. 
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spread over a bandwidth roughly twice that of the transponder output* 
These noncoherent products thus appear as noise within the Doppler filter 
bandwidth * 

One technique for generating the warning alarm is shown in Fig. 4, 

2 

utilizing a warning parameter designated as (R x)» In this technique, 
the Doppler signal is differentiated and detected to provide a voltage 
related to "time to collision." A threshold circuit is then used to give 
an alarm when this voltage exceeds a certain level . 


e_ = A sin go , t e 0 = K- Ago xos GOjt e 0 = K^Aoo 

1 d L Id d 3 2 i 



A is proportional to R 


-3 


a) . is proportional to closing velocity (V ) , 

CL 2 c 

hence voltage e^ is proportional to (R x)“^ 


Fig. 4. Technique for generation of alarm 
based on R^x criteria. 


Two of the systems have been constructed and installed in DC-4 
aircraft for flight test evaluation. 

In the systems subjected to flight tests, signal processing techniques 
were used that permitted derivation of voltages proportional to range, 
closing velocity, and rate of change of closing velocity, as well as 
combinations of these parameters. Both analog and digital data process- 
ing techniques have been evaluated during the flight test experiments . 

Additional details of the system operation are given in reference [1] 
and [3], and in the following sections of this report. 
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B, SYSTEM SIMULATIONS 


Several simulation programs have been developed, each somewhat 
similar, but differing in input data, output, and in certain internal 
details. A general description of each of these programs follows: 

1. Linear path simulation : 


inputs : 

positions and velocities of up to 
15 air crafty system parameters, and 
time intervals. 

outputs : 

geometrical parameters and warning 
criteria, transponder and receiver 
power levels, frequencies, alarm 
status, signal-to-noise ratios, 
saturation condition, interference 
levels and output voltages vs time. 
(For all transponders and one se- 
lected receiver). 

Atlanta data 

simulation: 

inputs : 

Atlanta radar data, system 
parameters, track selected for 
analysis, and time intervals. 

outputs : 

same as in (1) 

Flight test 

simulation: 

inputs : 

Radar data from Wallops FPQ-6 and 
FPS-16 radars, system parameters, and 
experimental flight test data. 
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outputs : 


Geometrical parameters and warning 
criteria* system outputs vs time 
as in (1), receiver voltages 
corresponding to geometrical 
parameters and warning criteria* break- 
down of system errors, and comparison 
of geometrical (ground radar) 
calculated (simulated warning system) 
and experimental (flight test data) 
measurements . 

4 . Saturation evaluation : 

inputs: Atlanta radar data and system 

parameters 

outputs: Statistical data on system power levels * 

alarm condition, and saturation 
condition, assuming full simulated 
warning systems on all aircraft in the 
data base a 

A general flow chart for the system simulation program is shown in Fig* 5* 
This chart indicates the basic program which is used, with slight modifica- 
tions, in the programs described above. 

C* ANALYTICAL MODELS 

In the phase I report [1] analytical models are documented that provides 
for calculations of signal, noise, and interference levels at various points 
in the transponders and receivers* Additional analytical and emperical 
models have been developed in phase II of the study to provide a more 
realistic representation of the behavior of the systems as presently con- 
structed* 
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CROSS-PRODUCT POWER MATRIX 
UNDESIRED SIGNAL POWER MATRIX 

SIGNAL POWER AND FREQUENCY FROM EACH TRANSPONDER 
SIGNAL-TO-NOISE RATIO MATRIX 

MEASURED RANGE, RANGE RATE AND WARNING CRITERIA 
ALARM STATUS 


Figure 5. Simulation flow diagram 
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The major refinements incorporated in the phase II simulations are 
summarized as: 

1. More accurate representation of the system antenna 

patterns using a Fourier analysis of the 
experimental patterns. 

2. Representation of the receiver amplitude-frequency 

response by empirical functions fitted to 
experimental calibration data taken at LRC. 

3. Development of a mathematical model of the output 

of a limiter-product detector circuit with 
several input signals 

4. Consideration of each transmitted signal separately 

instead of assuming equal path losses and 
antenna characteristics for both transmitted 
frequencies. 

5. Incorporation of time response in the receiver 

calculations in accordance with time constants 
used in the actual systems. 

6. Provision for the simulation of the analog data 

processing accomplished in the receiver based 
on experimental calibration data, 

7. Incorporation of an error breakdown giving calculated 

measurement errors as referenced to an ideal 
receiver and to experimental flight test data. 

Details of the above refinements are given in Appendix D (Receiver 
Data Processing), Appendix B (The Output of a Product-Detector Limiter 
Circuit) and Appendix G (Representation of Antenna Patterns). 
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IV. SATURATION AND INTERFERENCE STUDIES 


A. GENERAL 

The radar traffic data from the Atlanta terminal , as mentioned in 
the introduction and analyzed in detail in reference [2] , permit studies 
of particular aspects of the system performance by computer simulation. 

A detailed warning system can be simulated on each aircraft in the data 
base* and statistics determined on performance measures of interest. 

In the study of the saturation problem, the system parameters of 
interest include the signal, interference, and noise levels at various 
points in the transponder and receiver. A good measure of the severity 
of the saturation problem is given by a statistic such as "the average 
percentage of flying time in the terminal area that a certain power 
level exceeded a specified level." This statistic also provides an 
unbiased estimate of the probability that a randomly selected system will 
have a power level exceeding the specific level at any particular instant 
of time. 

Another statistic of interest is"the average percentage of flying 
time in the terminal area that a certain power level exceeded a specified 
level while the system was in an alarm status ." This statistic indicates 
the severity of both the multiple hazard and interference problems. For 
example, if the second largest signal power level at a given receiver 
exceeded the alarm threshold of the receiver for a large percentage of 
the time, a multiple target problem would exist. As another example, if 
the signal-to-interf erence ratio for the largest signal at a given 
receiver (while in an alarm status) is small for a large percentage of 
the time, an interference problem exists. 

Thus, in the following, results are presented indicating the above 
mentioned statistics on various power levels for hour 11 of the Atlanta 
data base. During this hour, 68 aircraft were present with a total flying 
time of 13.6 hours [2]. 
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Bo TRANSPONDER CHARACTERISTICS 


1. Mixer Saturation 

The linear turn-around transponder deviates from linearity due 
to two major effects* One effect is that of exceeding the dynamic range 
of the balanced mixer circuit that effectively multiplies the two 
incoming signals and derives the low-level output signal (see figure 2)« 
The other effect is the deviation from linearity of the output linear 
amplifier chain* Experimental data indicate that with one set of input 
signals j the balanced mixer is the component that tends to saturate first 
as the input power levels are increased* 

A theoretical study of mixer saturation (see reference [1]) indicated 
that the mixer saturation was a function of the total power input to the 
mixer s and that the effect of transponder saturation could be represented 
satisfactorily by a models 
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Pb 


ii 


1 + K 


N 

l 

i=l 


(Pa 


ij 


P V 


( 1 ) 


where 


Ps 


ij 



Pb 


ij 


F 

x 

K 


power output of transponder j due to transmitter in 
aircraft i 

power input from aircraft i at transponder j at 
frequency a 

power input from aircraft i at transponder j at 
frequency b 

overall gain constant for transponder 

experimentally determined constant approximately 
equal to the gain of the transponder from the 
input to the mixer, divided by the mixer bias 
power (units of MW P ) 
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The denominator of eq. (1) has been designated as the "saturation factor" 
for the transponder * Notice that this term is also the nominal gain 
divided by the actual gain for a given signal, or 


Sat, factor (S v ) 


nominal gain 
actual gain 


( 2 ) 


Figure 6 shows the results of the simulation study indicating the 
average percentage of time that an aircraft transponder "saturation 
factor" exceeded a level specified along the horizontal axis. For example , 
a value of 1.3 was exceeded only .6% of the time; this corresponds to a 
gain reduction of only 1.1 db. 

Inspection of the curve indicates that the median gain reduction 
factor is 1.05 and that there is a low (= .01) probability of exceeding 
1.4. These levels are acceptable, hence transponder mixer saturation 
does not appear to be a problem under hour 11 (congested) flight 
conditions, and with the system parameters (gains, bandwidths, etc.) 
used in the construction and simulation of the system. 

2. Signal and Interference Levels 

Another factor of importance in judging the transponder performance 
is the level of the cross modulation products in the transponder resulting 
from multiple interrogations. If n aircraft are interrogating the 
transponder a total of n(n-l) cross products are generated. Because of the 
random modulation of the transmitted signals ? the cross products are 
spread over a bandwidth determined by the frequency deviation of the two 
oscillators in the transmitter (see spectral density in Fig. 3). 

Curves indicating the average percent of time that the cross product 
power and the noise plus the cross product power exceeded a level greater 
than a specified value are shown in Fig. 7. The power levels are referred 
to the transponder output (30 dbm maximum output). As may be seen from the 
plots, the noise plus cross product power level referred to the output 
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PERCENT OF TIME 



TRANSPONDER OUTPUT POWER (P k ) , DBM 


Figure 7. Average percent of time that the power output of the 
transponder is greater than for various outputs. 


17 



rarely exceeds 9 dbm. The median value of noise plus cross product power 
is approximately -6 dbm. 

Figure 7 also plots curves indicating the average percent of time 
that the total output, the largest signal output, and the second largest 
signal output exceeded the value specified on the horizontal axis. Note 
that the total output power is well below the 30 dbm maximum value for all 
cases. 

From the curves of Fig. 7, we conclude that the transponder is 
operating well within its design values under hour 11 conditions. For 
conditions existing during this hour the cross product power is not an 
excessive percentage of the total output power, and the probability of 
exceeding the design maximum output of the transponder is very low. 

C. RECEIVER 

1. Signals in Receiver 60 MHz IF Amplifier 

Figure 8 plots the largest, second largest, etc. up to the fifth 
largest signal in the receiver IF amplifier. Again these curves represent 
the percent of time that the power level exceeded the value along the 
horizontal axis, for hour 11 data. These power levels are referred to the 
input of the IF amplifier (-85 dbm corresponds to a target at a range of 
5 miles). 

The statistical analyses of the same hour of data reported in the Phase 
II Interim Report [2] indicated that 60 percent of the time there would be 
at least one aircraft within 5 miles of a randomly selected aircraft. The 
data in Fig. 8 indicate that the signal level from the closest aircraft 
exceeds the power level corresponding to 5 nautical miles only 30 percent 
of the time. The difference between these percentage values is evidently 
due to the filtering by the system antenna characteristics. 

The receiver noise level is approximately -104 dbm* hence the signal- 
to-noise ratios are high for a large percentage of the time. The largest 
signal exceeds the noise level 75 percent of the time. The largest 
of the curves indicate that a dynamic range in the receiver IF of 45 d£> 
should be adequate to prevent saturation, with high probability. 
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2. Receiver Doppler Filter Output 

The Doppler filter in the receiver is assumed to have a 3 KHz 
bandwidth and an ideal frequency response from 100 to 3,000 Hz. 

Figure 9 plots the percent of time that various power levels at the 
receiver Doppler filter exceeded the level plotted along the horizontal 
axis. The power levels are referenced to a level of -76.4 dbm, 
corresponding to a target at 4 nautical miles and closing at a velocity 
of 240 knots. This reference level was selected so that a 1,000 cycle 
Doppler signal from a target at 5 miles would provide a -85 dbm power 
level that was used in the preceding plot (Figure 8). Even though 
the signals are referenced to a level corresponding to the IF amplifier 
input, it should be noted that the signals have been passed through a 
differentiating amplifier, hence the power levels are a function of the 
Doppler frequency . 

Figure 9 also plots "average percent of time" curves for the cross 
product power levels generated in the transponder and the "undesired 
signal power". This latter power represents signals that arise from 
transponder signals that are returns to other receivers in the population 
of aircraft. The cross product power and undesired signal power both 
represent noise like signals spread across the doppler bandwidth. These 
power levels have been combined with the noise power level in a composite 
curve designated as the noise plus interference power in Fig. 9. 

We note from the curves that there is a low probability that the noise 
plus the interference power will exceed the threshold level of the 
detector. Thus, alarms caused by interfering signals should be negligible 
under the conditions represented during this hour of the data base. The 
cross product power term is the most signficant interference source. By 
extrapolation of the noise plus interference power level curve, we can 
estimate the probability of a false alarm due to interference as 4 x 10 

3. Receiver Signal Levels Under Alarm Conditions 

The preceding curves on the receiver power levels have considered 
the overall receiver operation during flights in the terminal area. 

A major interest, however, is the signal level during an actual alarm 
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condition. That is, we are interested in the noise and interfering 
signal levels while the system is in an alarm status. 

Figure 10 plots the percent of time that the received signal 
exceeded a level R^. for the receiver in an alarm status. For these 
curves, a hazardous situation (alarm condition) is defined by 


R^t < 1000 nm^-sec (3) 

where R is the relative range and t is the approximate time to closest 
approach. The received power level corresponding to the above geometrical 
definition is -76.4 dbm. In terms of received power levels and doppler 
frequency, the expression for the differentiated power is given by equation 
(4) 

P d - -79 + 20 log dbm (4) 

where is the differentiated power and f^ is the doppler frequency in Hz. 

As may be seen in Figure 10, the second largest signal and third largest 
signal while in an alarm status are less than the signal causing the alarm 
with high probability. The system under hour 11 conditions indicated an 
alarm 2.4 percent of the time while the second largest signal exceeded the 
threshold level approximately .04 percent of the time. 

Figure 11 indicates the percent of time that the receiver signal to 
noise plus interference power ratio is greater than the level plotted along 
the horizontal axis for threat situations only. As may be seen from the 
curve, the probability that the largest signal to noise ratio exceeds 20 db 
is on the order of .9, given that the receiver is in an alarm status. 

No cases were observed in which the signal to noise ratio was less than 
10 db for the signal causing the alarm. 

From Figs. 10 and 11 we can conclude that, when in an alarm statuq with 
the threshold condition defined by eqs. 3 and 4, the interfering signals and 
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Figure 11. Average percentage of time that the receiver signal to 
noise and interference power ratio is greater than SNR^. 
for alarm situations only. 
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noise will not be at a level such that they interfere with the proper 
operation of the warning system. The signal to noise plus interference 
ratios during an alarm are sufficiently high such that additional informa- 
tion, such as range acceleration, can be derived with the signal processing 
equipment . 

4. Typical Signal Spectra in While in Alarm Status 

It is informative to note in detail the power levels and 
frequencies of the five largest signals (with closing doppler) from which 
the statistics of Fig. 10 and 11 are obtained. Table 1 shows a partial 
list of these five largest signals during alarm conditions. The Table 
indicates the power levels, signal to noise plus interference ratios, and 
doppler frequencies as shown in the Table heading. The specific cases shown 
in Table 1 were selected at random from the total of 131 alarm cases 
observed in the hour 11 data base. 

It should be noted that in most cases the fifth largest signal is 
considerably below the receiver noise level hence, the data indicate 
approximately the situation that would be observed on a audio signal 
analyzer looking at the differentiated signal from the doppler filter in 
the receiver, (i.e. usually less than 5 signals are observable) 

D. RECEIVED SIGNAL AND INTERFERENCE LEVELS VS TIME 

While the operation of the system is adequately explained by the 
statistical curves given in the preceding sections, it. is interesting to 
examine the time history of specific receiver operation. For this reason, 
a particular receiver was selected from the "threat selection tables" 
given in reference 1. This track was selected because a situation occurred 
in which an alarm should have been noted on the receiver. The track 
selected was track 98 in the hour 11 data base, which had a total flight 
time of approximately 11 minutes and was arriving at the terminal. 

Figure 12 plots the received power levels versus time for the largest 
signal received, the largest differentiated signal, and the noise plus 
interference power level versus time. For this particular simulation run. 
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Figure 12. Signals and noise power received at receiver on aircraft 
of track 98 (hour 11) . Time base is the time from the 
start of the track. 



SIGNALS REDUCED IN THIS INTERVAL BECAUSE 
A/C OF TRACK 98 WAS HEADED AWAY FROM ALL 
A/C IN THE DATA BASE (SEE FLIGHT PATH IN 
FIGURE 13.) 
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Figure 12. (Continued) 




a minus 67 dbm power level represented a threshold corresponding to 
2 2 

R t « 1,000 nm - sec (note that the threshold level is different from that 
used on the statistical curves, however the geometrical warning criteria 
is the same as given in eq. 3)* Note that the noise plus interference 
power level is well down from the signal levels throughout the whole track. 
Three separate alarm periods are indicated by the shaded area where the 
largest differentiated signal exceeds the threshold level. 

At 6 1/2 minutes after the start of the track the signals appear to 
drop out and all power levels reduced to the receiver noise level. The 
reason for this dropout was examined in detail and determined to be due to 
the fact that the aircraft of track 98 was heading away from all other 
aircraft in the data base during the time interval of the dropout. Figure 
13 plots the flight path of track 98 so that the position on the track can 

be correlated with the signals received. At the time of the signal dropouts 

the aircraft track 98 was headed away from the airport at a range of 
approximately 12 miles. 

Also indicated on Fig. 13 are the tracks of aircraft that caused the 
receiver of track 97 to indicate an alarm. The first two alarms were 
received just after 4 minutes in the track when the aircraft observed other 
aircraft on the final approach to the runway (track 87). The third alarm 

occurred due to a crossover of track 97 just after 6 minutes into the track. 

At the time the alarm was received on track 98, the approximate time to 
closest approach was 33 seconds and the projected miss distance was .75 
n. mi. Track 98 was at 2700 ft altitude and descending, while track 97 was 
at 1700 ft. and flying level. At the point of maximum signal, the slant 
range was .81 n, mi. and the approximate time to closest approach was 24 
seconds. The altitude separation at the point of closest approach was 
somewhat less than 250 ft. 

From examination of this specific case, which can probably be considered 
typical, we see again that the interference levels are not severe and that 
the differentiated signal does give a good indication of a potentially 
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Figure 13. Flight path of track 98 (hour 11) indicating alarms received. Numbers in 
parentheses are times in minutes from the start of the track. 




hazardous situation. The first two alarms received at the receiver occurred 
during a turning maneuver by the tracked aircraft and hence did not 
represent a hazardous situation. The last alarm, however, would certainly be an 
encounter of the type that would be considered potentially hazardous. 

E. DISTRIBUTION OF POWER LEVELS FOR IDEALIZED SYSTEMS 

For system design, it is important to be able to estimate the range 

of signal levels to be received from the desired target as well as that 

from other targets which may cause interference. While this has been 

accomplished by direct simulation of the specific LRC system using the 

data base, certain general results from the statistical analyses of 

reference [2] prove useful for estimation of dynamic ranges and interference 

levels for other systems under consideration. 

A histogram of the distribution of the range to the closest aircraft 

from hour 11 data is hown in Fig. 14. This histogram represents an 

approximation to the probability density function of the probability of 

observing, from a randomly selected aircraft, another aircraft within a 

range increment AR. The histogram data, when plotted on lognormal 

probability paper, indicates that the range distribution over the region 

0-10 n. mi. can be approximated by a lognormal density function with mean 

2 

(y^) and variance (a^ ) of: 


y ^ = In 4.0 

(5) 



where the basic units are nautical miles. 

A plot of the lognormal distribution corresponding to these parameters 
is shown on Fig. 14 for comparison with the histogram data. Thus, we have. 


f x (R) 


Exp {- -t- (In R-p ) 2 } 

/2tT c^R 2 o ± 


( 6 ) 
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where f^(R) is the probability density function of the range to the kth 
closest aircraft from a randomly selected aircraft under hour 11 conditions. 

Similarly, for the second closest aircraft to a randomly selected 
aircraft, we have the distribution shown in Fig. 15. The mean and 
standard deviation for these data are 


p 2 = i n 6.0 

°2 ‘ - 44 


(7) 


2 

and f ? (R) ~ A (In 6.0 , .44)* where A (p*a ) designates the lognormal distribu- 

Z 2 
tion with mean jj and variance a » 

For the third* fourth* and fifth closest aircraft* we have from the 
data* 



f 3 (R) = 

A (In 8.0, .38), 

(8) 


f 4 ( R ) = 

A (In 9.6, .25), 

(9) 

and 

f 5 (R) = 

A (In 11.4, .22). 

(10) 


The fits of the lognormal functions to the data are shown in the plot 
on lognormal paper* Fig. 16- As may be seen* the fits are not perfect* 
however* no other analytical distribution was found that provided a better 
fit* The distributions investigated included the Rayleigh* Rice* Poisson* 
and normal distributions. 

Some theoretical justification can be given for the lognormal 
distribution in that it arises under conditions when a change in a 
variate is a random proportion of the previous value of the variate [5]. 
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Figure 16. Relative range distributions plotted on lognormal 
scale to indicate fit of data to the lognormal 





Thus s if we assume that pilots fly such that the relative closing velocity 
is proportional to the relative range, with a random proportionality 
constant ; we have 


AR. 



(ID 


where the j subscript indicates the j_th time instant and At is the time 
increment*. Rearranging Equation (11) and summing over j gives 


n 


l 

j=i 


AR. n 

— L = y e 

At R. £ 

3 j=l J 


( 12 ) 


If each increment is small, and taking At = 1 without loss of generality. 


1 n AR. 

“ 3-1 ^ : 


R 

• n 

dR _ - 

— = In R - In R 
R n o 


(13) 


where R q is the initial range and R^ is the range at j = n. Thus, 

In R = In R + e- + e 0 + .... e (14) 

n o 1 2 n 7 

By the central limit theorem, In R^ will be asymptotically normal (since 
it is the sum of a large number of small random effects), hence R^ will 
be distributed lognormally. Note that the same result is obtained if we 
assume that the change in relative range is randomly proportional to 
the existing range. 

The lognormal distribution has transformation properties that are 
extremely useful. These properties are summarized in a theorem from 
Aitchison and Brown [ 5 ]: 

2 3 , 

If X is A(p,a ) and b and c are constants (c = e ), 
then c x b is A (a + by, b^c^)." 
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Using this theorem* we can determine the probability density of the 
power received at a randomly selected aircraft i. To relate power 
received to the relative range, we assume isotropic antennas for both 
transmitting and receiving, and use the radar range equations: 


-6 V 

R case: P. = -z — 

J (4tt) RL 


b , 6 

- Kj /» 


, P t X 2 a. . 

R 4 case: P. = — Y — 1 * L “ “ K r /R 4 

J (4tt ) 3 R 4 L * 


where 


j 

r jb 


= power received at aircraft i from aircraft j 
= power transmitted at aircraft j 
= transmitted wavelength 
~ system losses 

= radar cross section of target j 
= power transmitted at i 

-6 -4 

= superscripts indicating R and R cases • 


b 

Thus, if f^XP) denotes the probability density function of the power 
received at i from the kth closest aircraft to aircraft i, we can write, 
using the theorem above and the previously found distributions; 


R-6 case: 


f^(P) = A (In kJ - 6 In 4.0, 9) 
fk(P) = A (In - 6 In 6.0, 15.8) 


(15) 

(16) 


(17) 

(18) 
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f^(P) = A (In - 6 In 8.0, 13.7) 


(19) 


f£(P) = A (In kJ - 6 In 9.6, 9) (20) 

f^(P) = A (In Kj - 6 In 11.4, 7.9) (21) 

-4 

Similar results are easily found for the R (two-way path) case if desired. 

Care should be taken in using these density functions since the 
f^(P) for K > 1 are conditional density functions, with the condition being 
that there are K-l aircraft closer to aircraft i. For isotropic antennas, 
and if all are equal, this is equivalent to saying that K-l aircraft 
provide a larger power level at aircraft i. 

A further useful property of the lognormal distribution is that: 

A(x) = N(ln x) (22) 

where N indicates a normal distribution. Using this property and the 
previously mentioned transformation theorem, we can write 

10 log x = 4.34 In x 

therefore, A(P^*^) = N(ln P^’^) = N(1Q log P) (23) 

Hence, the power in db or dbm is normally distributed . We apply the 
4 34 b 

transformation P * to the f R (P) functions using the transformation 
theorem as follows: 


if 

f£< p ) - A(p k , c£) 


then 

f k (p4 ‘ 34) " A < 4 - 34 v 18 - 8 ° k 2) 


and 

fJ(P db > = N(4.34y k , 18.8o k 2 ) 

(24) 

where P^ = 10 log P. 

-4 

Note that if the R (two way path) 

range equation is 


used, the power in db or dbm is still distributed normally. 
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note 


In order to illustrate the results, we assume all are equal, and 

that represents the power received at aircraft i from another 

another aircraft at a range of one nautical mile* 10 log is this 

power in db or dbm depending on the reference power units* If we 

b 

designate this power level as t * ie expressions for the density 


functions become: 


fb l< p db> 

- B ' S db - 

- 36.0, 

169) 

(25) 

£ 2 <P db' ) 

- S<S db ' 

- 46.8, 

298) 

(26) 

£ 3 (I W 

‘ B(S db - 

- 54.3, 

257) 

(27) 

f 4 (P db> 

- K(S db ' 

- 58.8, 

169) 

(28) 

£ 5 (P db> 

’ »< S db ' 

- 68.7, 

148) 

(29) 


The probability that the power received from the kth closest aircraft 
will exceed a level £ is determined from 


Prob* (P 


db 


0-1 


W dP db- 


(30) 


Curves indicating this probability for K = 1, . .., 5 are shown in 
Fig, 17. 
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(31) 

(32) 

(33) 

(34) 


f 5 (P db ) = N(S db - 45.8, 66.0) (35) 

T 

where S ^ is the reference power level (i.e. power received from target at 
one n. mi. ) * 
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V. FLIGHT TEST EVALUATION OF RANGE MEASUREMENT ERROR 


A, GENERAL DESCRIPTION 

Two of the warning systems constructed by LRC personnel were installed 
in a pair of DC-4 aircraft and subjected to flight tests* The major 
objectives of the flight test series were to evaluate measurement errors, 
to determine the sources and magnitudes of errors, and to determine the 
system performance under actual flight conditions* 

During the flight test operations at Wallops Island, Virginia, the 
positions of the two aircraft were measured by the Wallops FPQ-6 and FPS-16 
radar systems* The radar data obtained during the flight test was then 
processed by Wallops personnel to derive relative data between the pair 
of aircraft. This data was furnished to RTI on digital magnetic tape to 
provide an input to the flight test simulations and to permit comparison 
of simulated, experimental, and geometrical measurements. 

The flight paths flown were a series of converging pairs with various 
predetermined miss distances, both horizontally and in altitude. In all 
cases, the aircraft were separated in altitude to prevent any danger to the 
aircraft. Table 2 summarizes the flight test designations, the altitude 
separation, the attempted and actual miss distances, and the magnetic 
headings of the aircrafts as provided by the pilots* The aircraft air 
speeds were on the order of 140-160 knots for all flights, and in all cases 
straight line flight paths were attempted by the pilots. 

Two computer programs were developed to process the flight test data. 

One program designated as GE01 processed the Wallops radar data and generated 
the aircraft flight paths along with geometrical parameters and values 
of warning criteria calculated from the geometry* The second program was 
the modification of the simulation program as described in Section III* 

This latter program is discussed in detail in the following section. 

B* FLIGHT TEST SIMULATIONS 

During flight test operations, the relative azimuth angle, relative 
elevation angle and relative range between the two aircraft are measured 
with the ground radar system* Calculations made using the ground radar 
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Table 2. Flight test parameters . 
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measurements are designated as geometrical values. Also, during flight 
test operations the received signal voltages are recorded to give the 
values of parameters as measured with the warning system. Calculations 
made using these latter measurements are designated as experimental 
values. Parameter values calculated using the simulation program are 
designated as calculated values. 

Figure 18 shows an example of the computer printouts obtained at 
5 second intervals throughout each of the flight test. The first line on 
the printout provides the flight identification, the elapsed time since 
the start of the flight, and the time of day (gmt) . The next section on 
the printout indicates geometrical parameters calculated directly from the 
ground radar data. These parameters include: 

(1) xyz coordinates for both aircraft 

(2) the relative range between aircraft 

(3) closing velocity in knots (V ) 

(4) the rate of change of closing velocity (V ^) 

• <* * 

(5) the normalized range acceleration, R/R (v ) 

(6) the normal velocity component (V ) 

(7) the time to closest approach (T) 

(8) the approximate time to closest approach (x) 

(9) the approximate miss distance (D) 

(10) the exact miss distance (R ) 

2 ° 

(11) geometrical value of R x (g) 

(12) the value of modified tau (x ) 

(13) the relative azimuth and elevation angles for 

both aircraft 

The next section on the printout sheet gives the transponder outputs for 
system 1 and system 2. Following this list the gain variations of each of 
the system antennas are listed. The gain variations are referenced to a 
0 db level (head-on conditions) . 
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Figure 18. Example of computer printout for flight test data. These printouts 
were obtained at five second intervals during each test. 



The receiver outputs are listed for both systems* These outputs 
include signal levels, signal to noise plus interference ratio, AGC 
voltage, doppler frequency, and the alarm status of the receiver* 

The next listing on the computer printout list geometrical, experimental, 
and calculated values of various geometrical parameters and warning 
criteria* The calculated values are given in geometrical units and also 
as voltages measured at the various receiver outputs* Also, in this 
section the percent errors are calculated for the difference between 
experimental and geometrical values, calculated and geometrical values, 
and experimental and calculated values, A description of the expressions 
used to calculate the values in this section is given in Appendix D . 

The last section on the computer printout gives an error breakdown 
for 3 parameters measured by the receiver; range, tau, and beta* The first 
column in this section repeats the percent error between the experimental 
and geometrical (ground radar) values of the parameter* The second column 
indicates the error due to amplitude variations in the receiver doppler 
filter* The third column gives the percent error due to deviations of the 
detector characteristics from linearity* The fourth column indicates the 
error and measurement due to gain variations of the elevation patterns 
(transmitter, transponder, and receiver). The fifth column indicates the 
percent errors due to the complete antenna patterns, again considering all 
patterns in the system* The sixth column is a fixed error source due to 
gain calibration errors in the flight test system* That is, the calculated 
or simulated values of range are based on exact gain calibrations that 
provide a -85 dbm signal level at a range of 5 miles* During flight test 
it was determined that the overall gain of the warning system was down 
somewhat, and this error source reflects the value of range and analogous 
to the deficiency in system gain of the flight systems* 

The last two columns in the error breakdown section indicate residual 
errors between the experimental and geometrical values that are unaccounted 
for by the error sources previously considered* The first residual 
column does not take into account variations in the azimuth patterns of 
the antennas, while the last residual column is the percent error remaining 
after effects of the doppler filter, detector slope, and the total 
antenna patterns are removed* The percent error remaining in this last 


45 



column 3 then* is due to system errors not subtracted out in the calcula- 
tion (e*g 6 saturation of the transponder) and factors not considered in 
the simulation such as aircraft deviations from straight line, level 
flight. 

It should be noted that the systems were calibrated over a range from 
1 to 5 nautical miles, hence relatively large errors outside of this range 
were expected and are felt to be unimportant in system evaluation, 

C* FLIGHT TEST RESULTS 

Extensive analyses of the flight test data have been conducted by LRC 
personnel, using the results of the flight test simulation program to 
determine contribution to the measurement errors of various system components. 
Thus, in this section the major emphasis will be on documentation of the 
flight paths and comparison of the simulated and experimental relative range 
data. Since the relative velocity between aircraft is simply a frequency 
measurement on the doppler signal, errors in measurement of this parameter 
were small. 

For each flight test as listed in Table 2, the following plots have been 
made : 

(1) the flight path of the aircrafts. 

(2) the geometrical, calculated and experimental 
range measurements vs. flight time. 

(3) the percent difference between calculated 
and experimental range measurements plotted 
vs. the value of geometrical range. 

These curves indicate the accuracy of the simulation of the actual flights. 

It should be recalled that the residual differences between the calculated 
and experimentally measured values of range are most likely due to variations 
of the aircraft from straight line level flight. Computer printouts listing 
the data at the approximate time the alarm was received in the simulated 
system are given in Appendix F, 
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Inspection of the plots indicates that for most flights the difference 
between simulated and experimental values of range was within 15 percent 
over the calibration range of the equipment . An exception to this agree^- 
ment was found in flight 8C 9 where the experimental and calculated values 
disagreed by as much as 80 percent. Although the cause for this 
discrepancy has not yet been determined 9 it is felt to be due to a mistake 
in calibration of the experimental data in the digitizing process® 
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VI, MISCELLANEOUS STUDIES 


A, GENERAL 

During the course of the contract, several studies were undertaken 
that have general applicability to the collision-warning problem. These 
studies include: (1) common definitions of warning times for various 

system warning criteria, (2) the possibility of alarm suppression due to 
relative range acceleration threshold logic, (3) comparison of warning 
times of various systems with equal alarm probability, and (4) warning 
times required for escape maneuvers. 

Other studies were concerned with extension of the "modified tau n 
warning criteria developed by Collins personnel. Reference [4], The 
extensions involved consideration of relative aircraft accelerations as a 
random variable and the introduction of additional constraints on aircraft 
acceleration vectors. Also, general results were derived for determination 
of the DC output of a product detector-limiter circuit with multiple 
input signals. Results of these latter studies are given in Appendices A 
and B respectively, while the general studies on system warning times are 
presented in the following. 

B. DEFINITION OF WARNING TIMES 

Since some system warning criteria and warning times are derived 
based on the assumption of nonaccelerating flight, while others are based 
on accelerating flight with acceleration constraints, it can be misleading 
to speak of a warning time attained without further elaboration as to 
the assumptions used. 

To provide a common basis for comparison of the protection afforded 
by the various systems, we define warning times as follows: 

t^ - minimum warning time for a nonaccelerating 
co-altitude collision with a maximum 
relative velocity of V . (i.e, this is the 

time-to-collision assuming a roll-out to 
linear flight at the instant the alarm is 
received, and a subsequent collision) 
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t k = minimum warning time achieved for a co-altitude 

threat under the worse possible conditions 

within the assumed velocity (V ) and accelera- 

m 

tion (U) constraints. (i.e. this is the time- 
to-collision assuming flight paths subsequent 
to the alarm at the maximum relative accelera- 
tion and/or the maximum relative velocity until 
collision occurs) 

V mv ~ warn ^ n S for a nonaccelerating collision 

occurring in the vertical plane with a maximum 

relative rate of descent or ascent of h • 

m 

In Section VI-C, inhibition of the alarm by range acceleration 
discrimination is considered. This inhibition will be shown to occur 
for sets (S^) of particular trajectories, hence to distinguish warning 
times for these cases, we define, for systems using R* measurements for 
d is cr imina t ion : 


t n st mh ” same as P rev ^ ous definitions, except trajectories 
are assumed not in S- 

d d 

t , t , = same definitions as t , t , , except that the times 
n mh n mh r 

are calculated using the worst-case member of the 
set of trajectories (i.e. for the trajectory 

that provides the smallest warning time.) 

Expressions for the times t^, t and t are easily derived. Let 
R(t ) and R(t Q ) be the relative range and closing velocity at which a 
given system gives an alarm. The time t^ is then 


t « min 
n 



over C 


( 36 ) 


where C is the contour in the R, R plane defined by the system warning 
criteria. 
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The minimum warning time for a system that provides an alarm at 
range R(t Q ) and closing velocity R(t Q ), assuming subsequent relative 
accelerations no greater than U is derived in [4] (i,e. the modified 
tau expression) 

-R(t ) + [R 2 (t o ) + 2UR(t Q ) ] 1/2 
t l = U 


( 37 ) 


where t^ designates this warning time. For the case where the trajectory 

providing the alarm time t^ does not exceed the velocity constraint, we 

have for the time t , 

m 


t ffl = min { t 1 [R(t Q ), R(t Q )] over C (38) 

where C is again the contour in the R, R plane defined by the system warning 
criteria. If the trajectory exceeds the velocity constraint, (i.e. if 
R(t Q ) - UT 1 > V ) graphical or approximation techniques can be used to 
determine t^, assisted by the phase plane plot of trajectory slope isoclines 
given in Appendix C. 

For those systems using altitude discrimination of the form | AA | < AA^, 
the minimum warning time for closing paths in the vertical plane, assuming 
no acceleration in this plane, is 


t 

v 



(39) 


Thus, by assuming a maximum closing velocity V^, a maximum rate of 
ascent (or descent) of h^, and a maximum relative acceleration U, the 
defined warning times achieved by the various systems may be determined 
to provide a common basis of comparison of the protection afforded by the 
systems. 
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Table 3* Warning times t n and as defined in Section VI-B for various system 

types. The maximum closing velocity assumed is Y m * the maximum relative 
acceleration is U and the maximum relative rate of descent or ascent is 
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Notes If altitude discrimination (|AA| < e AA^) is used* the minimum warning time is the value 
from the expression listed^ or AA^/h^, whichever is smaller. The minimum warning times 
may also be determined by R or normal velocity discrimination as discussed in Section YI-C, 



For the R x and tau systems, the minimum warning time approaches 
zero as the relative closing velocity approaches zero. However, a 
supplementary (over-riding) range threshold can be used to prevent this 
situation. Hence, in the work that follows, we will consider these 
systems with a supplementary range threshold of magnitude 

In Table 3 , expressions are given for the defined warning times for 
the various systems that have been considered, in terms of the warning 
thresholds and the assumed constraints V , h and U. The effect of range 
acceleration discrimination on the defined warning times is considered 
in the following section, 

C. ALARM SUPPRESSION DUE TO RELATIVE RANGE 
ACCELERATION THRESHOLDS 

For systems using discrimination based on R measurements (e.g. 
approximate miss distance thresholds assuming nonaccelerating flight) , 
the measurements will, of course, be erroneous if accelerating flight 
occurs. Under accelerating conditions, we must consider the possibility 
that accelerations due to maneuvers will effect the measurement in the 
worse possible way. These worse case conditions are considered for 
various systems in the following. 

Tau, Range, Normal Velocity System . In this system, normal velocity 
is estimated in accordance with [2] 



V = 1 77 

n R R 


(40) 


which was derived based on nonaccelerating flight paths. The system warning 
criteria are 


R/K < T k and 


< V 


nk 


( 41 ) 


or R < 




where the k subscripts indicate selected threshold constants. 
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For closings nonaccelerating flight* R is positive . For closing 
accelerating flight* however* R can be either positive or negative* 

If the magnitude of R is sensed, such that the threshold condition 
is 



V i and R/R < t. 
nk k 


( 42 ) 


and if relative accelerations no greater than U are allowed, the pos- 

i ** i 

sibility exists that measurements of |R| will equal U. If this is the 
case, an alarm will be suppressed so long as 


or until a range less than 


> V 


nk 



(43) 


(44) 


is reached. If this range is greater than the supplementary range 
threshold R^, and if the worse case trajectory (within the constraints) 
is assumed, the alarm would occur at the maximum closing velocity V . 
The warning time at this alarm point is determined by division of R^ 
by or, 



V 


t = 
n 


nk 

U V 9 
m 


V 


if 


nk 

U 


V 


(45) 


d d 

where t - and t were defined in Section VI. B. The situation is shown 
mh n 

in the phase space sketch of Fig. 34. 
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CLOSING VELOCITY, KTS 








As an example of the numerical value of warning time suppression, 
we have for V nk = 120 KTS, U - 1/2 g, and V m = 400 KTS, 


= .42 n. mi. 


( 46 ) 


Hence, unless the supplementary range threshold (R^) were greater than 
this value, we could collide 3.8 seconds after the alarm. 

If the sign of R is sensed, such that the threshold condition is 


R/R < t, and 
k — — 



< V . if R > 0 
nk 


(closing velocity decreasing) 

(47) 

R/R < t, if R < 0 (closing velocity increasing) 

K. 

or R < 


the situation is not quite so severe 9 although the alarm can still be 
drastically suppressed* We could still suppress the alarm so long as the 
relative trajectory is such that 


V nk < r (48) 

R — 

In the worse case, the path is such that the equality holds, and such that 
the trajectory passes through a point in the R, R plane; (T kV V- 
Since we allow acceleration magnitudes no greater than U, the path can 
suddenly change to one with acceleration of -U, and a collision can occur. 

The alarm will be given at the point of deviation from the path. Thus we 
are led to determine the point that provides the minimum warning time. 
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For a path defined by 


a © 

R = 



( 49 ) 


and passing through the point ( T k v m > v m )> the points in the R, *R plane 
defining the trajectory can be found by dividing Equation by R 
and integrating. 


• * V 

R _ nk _ dR 
; R R " dR 


Separation of variables and integration gives 


(50) 


•2 0 „ 2 R *2 

R = 2V In — + R , 
nk R. i* 

l 


(51) 


where R^, R^ are the initial values of range and closing velocity. 

The point of maximum hazard is given approximately by the point of 
minimum time to closest approach, or 


min {t} where x = R/R 


(52) 


(For accelerating flight, we should actually minimize t^ as given by 
Equation 37 , however, this minimization is more complex and graphical 
considerations indicate that the approximation of Equation 52 does not 
lead to significant error in finding the point of maximum hazard). The 
point on the path for maximum hazard is thus found as 


R = R. EXP ■=■ 
1 i 2 




( 53 ) 


/ 
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e 

where the initial conditions are R. = x, V and R. = V , for the worse-case 

x km x m 

trajectory. 

At the point (R^, R^) , the minimum time to collision is found from 
Equation 37 as 


d -*i + + 1/2 

fc mh U 


(54) 


which assumes an immediate trajectory change to increasing closing velocity. 

The situation is shown on the phase plane sketch of Figure 35 „ For the 

time t^, we have 
n 


t - R- /R- . 
n 11 


(55) 


If R^ is less than the supplementary range threshold R^, then the 
alarm would be given at the point 


R 






2 • 2 
2V. In -= + R. 
nk R. x 


] 


1/2 


(56) 


* 

j R|) in Equations 54 and 55 to 

determine the minimum times remaining until collision. 

A complication may arise in that the minimum time trajectory used in 

calculating the time given by Equation 54 may result in closing 

velocities exceeding the assumed maximum (V ) . For velocities on the 

m 

order of 2 miles * however, the actual time difference in suing Equation 

* 

versus the time for a nonaccelerating path (R^/R^) is small (e.g. 

~ 2.5 sec. for U = 1/2 g). 


Hence, if R^ < R^, we use (R^ 
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CLOSING VELOCITY, KTS 


400 



30 




R x, Range, Normal Velocity System , 
conditions are 


In this system, the threshold 


2 •• • 

R x <3, and R/R < v, 
k k 

(57) 

or R < R, 

and we will assume that only positive values of R are used. The alarm 
will be suppressed for trajectories such that 


R > y 1 R. 08) 

— k 

Values of range and closing velocity for the path such that the equality 
holds are found by integration as, 

i - R ± - Y k (R - R ± > > (59) 

where R. and R. are the initial points. 

IX 

In this case, the hazard increases continually for the worse case 
trajectory, and the alarm will be given when the supplementary range 
threshold is reached, or at, 


R i = R k 

(60) 

h = "k - V + V 

1/3 

The initial conditions for the worse-case trajectory are ( (g, V ) , V ) 

km m 

as shown in the sketch of Figure 36 . 
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CLOSING VELOCITY, KTS 


400 KTS' 


^POSSIBLE PATHS 
TO COLLISION ^ 


PATH FOR 
R/R < Yv 


V?OINT AT WHICH ALARM IS GIVEN 
FOR WORSE-CASE TRAJECTORY 

( V V 


•RANGE 


1000 NM -SEC 


RANGE, NAUT. MI. 


Figure 36. Phase plane plot for the R x-range-normal 
velocity system indicating possible alarm 
suppression due to a R/R < y^ threshold. 
For the sketch, y^ = .004 sec“l, and 
r2 t « 1000 nm2-sec are used. 



Modified Tau, Normal Velocity System . For this system, let us 
assume that normal velocity measurements are made by 


V = 
n 



R > 0 


V = 0 R < 0 

n — 


instead of the relationship valid for accelerating paths, 


( 61 ) 


V 

n 



(62) 


where V is the relative velocity. 

If measurements are made in accordance with Equation 61 , the 
alarm could be suppressed just as described for the tau-rang e-normal 
velocity system. Thus, equations derived for this system (i.e. Equations 
54 , 55 , and 56 can be used to determine the minimum alarm time under 

worse case conditions. 

The initial conditions for the worse case trajectory are 


Ux 


R. = V x , + 
l m mk 


mk 



(63) 


where the nomenclature has been previously described. 

D. WARNING TIMES OF VARIOUS SYSTEMS WITH 
EQUAL ALARM PROBABILITY 

In ref. [2], several systems that have been discussed in the literature 
were compared to determine alarm statistics that can be expected for 
approximately equal levels of protection (i.e. coverage in the R, R plane). 
We found unreasonably high alarm probabilities in some cases. In this 
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section* the type and range of measurements required and level of protection 

achieved are determined under constraints on the alarm statistics. In 

particular* the idealized systems simulated are compared using the 

statistical results in ref. [2] to determine warning times afforded for a 

given fraction of flying time in an alarm status. For the comparison, 

the results obtained from the most congested hour of data (hour 11) are used. 

For calculation of warning times, we use the results of Section VI-B. 

Table 4 indicates the measurements required, detection ranges, and 

warning times achieved for various ideal systems that were in an alarm 

status 1% of the time. The applicable threshold values were taken from 

Figures 3 through 11 in ref. [2]. The maximum detection ranges are based 

on a 400 KT closing velocity, as is the warning time (t ) as discussed in 

^ d 

the preceding section. The smallest of the warning times t t , and t , 

mh v mn 

indicates the minimum warning assured for closing velocities from 0 to 
400 KTS , for relative aircraft accelerations not exceeding 1/2 g, and for 
vertical rates of descent or ascent of 25 fps. 

2 

The zero minimum warning times in the t and R t systems have been 

eliminated in Table 4 by using a range threshold in conjunction with the 
2 

t or R t criteria. The supplementary range criteria was not investigated 

2 

directly (i.e. simulated with t and R t systems using the data base), 

however, we see from Fig. 3 that a 1 nm range threshold will not significantly 

increase the percent of time in an alarm status for those criteria without 

altitude discrimination. For systems using altitude discrimination of ± 500 

ft, a range threshold of 1.5 nm does not significantly increase the alarm 

probability. Hence, in Table 4 we have calculated a minimum warning time 
2 

for the x and R x criteria using a range threshold of 1 nm or 1.5 nm as 
applicable. The warning time calculations are made using the expressions 
listed in Table 3 and relationships in Section VI-C as applicable for the 
cases where R discrimination is used. 

E. WARNING TIMES REQUIRED FOR ESCAPE MANEUVERS 

A recent analysis of warning times (in the horizontal plane) required 
for a CAS has been conducted by Holt and Anderson [ 6 ]• In this analysis, 
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reasonable statistical distributions were assumed for altimeter errors , 
pilot reaction time, data processing time, and aircraft servo system 
delays • A computer simulation approach was then used to determine the 
statistical distribution of the time required for execution of an altitude 
separation maneuver to assure a 150 ft altitude separation. 

To execute the altitude evasive maneuver, where the pilot has been 
previously alerted and is in straight line, level flight, a warning time 
of 25 seconds is recommended as a result of the above study. 

For a roll-out alarm (i.e, if in a turn, stop turning) a warning 
time of 30 seconds is recommended. 

For a pilot alert, a warning time of 41 seconds is recommended. 

This allows an average pilot reaction time of 2,5 seconds and then 10 
seconds to level-off prior to execution of a possible evasive maneuver, 

(All times assume only one aircraft maneuvers). 

In the vertical plane, it is estimated by Perkinson [ 7], that for 

most aircraft, altitude rate can be reduced to zero within 15 seconds 

after warning. Hence, if descending or ascending aircraft level off after 
warning, a 21 second warning of an aircraft above or below would assure 

safe clearance of 150 ft for a maximum rate of descent (or ascent) of 

25 fps, assuming no vertical accelerations. 

For a pilot warning indicator that depends upon visual acquisition 
to evaluate the threat and to determine the required maneuver, further 
study is necessary to evaluate the time required to visually acquire the 
target and evaluate the threat. This is presently under study by the FAA * 

F. SUMMARY 

Comparison of the minimum assured warning time for systems that were 
found to be in an alarm status 1% of the time (Table 4 ) with the desired 
pilot alert time on the order of 40 seconds indicates that none of the 
systems provide this minimum warning time. Although the 1% alarm 
probability is quite arbitrary, it is felt to be a reasonable upper bound. 
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If we assume an average alarm duration of 10 seconds * Equation 
indicates that the probability of an alarm during a thirteen minute 
flight is ~ *5 under hour 11 conditions. Hence, we are reluctant 
to relax this constraint. 

The systems operating on range alone and range-range rate measure- 
ments only should be eliminated from further consideration as not feasible. 
The data clearly indicate that added discrimination in altitude and 
possibly in (true) normal velocity is desirable if the required warning 
times are to be achieved without excessive alarms under normal 
operational conditions. 

Discrimination based on range accleration measurements is effective 
in increasing the maximum warning time for fixed levels of alarm 
probability. Unfortunately, as indicated in Table 4 and Section VI-C, 
discrimination based on R measurements has an adverse effect on the 
minimum warning times achieved under certain worse-case trajectories, 

when accelerating flight is considered. The effect is particularly bad 

... 2 » , 
when measuring R/R in the R x-range system. Thus, R measurements do not 

appear desirable for use in alarm logic, although this measurement may 

prove useful for indications of the effectiveness of avoidance maneuvers 

after an alarm is received. 

Normal velocity discrimination is effective in the modified tau 
system, providing (under the conditions discussed) a minimum warning 
time of 32 seconds as compared to 26 seconds obtained without this form 
of discrimination. Derivation of true normal velocity, however, requires 
an exchange of velocity vector data between aircraft, and adds considerably 
to the system complexity. 

We thus conclude that, unless an unduly high alarm rate is accepted, 

the minimum measurements required are relative range, closing velocity, 

and altitude. With these measurements, minimum alarm times on the order of 

2 

30 seconds can be obtained using the modified tau or R x-range warning 
^criteria (see Table 4 ) under the acceleration and velocity constraints 
assumed. For longer warning times and/or reduced false alarms, it appears 
necessary to exchange velocity vector data between aircraft. 
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Under the velocity and acceleration constraints assumed, the 
required detection range is on the order of 5 miles for the modified 
tau - altitude system providing a 32 second warning time. 
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VII . CONCLUSIONS 


The simulation of the system using the most congested hour (hour 11) 
of the Atlanta data base indicated that interference levels from all 
sources were small in comparison to the signal from a potentially hazardous 
aircraft. Transponder saturation by multiple aircraft interrogations 
does not appear to be any problem with the system parameters now used. 

The probability, under hour 11 conditions, that the transponder saturation 
factor will exceed a value of 1.4 (i.e. a gain reduction of 1.45 db) 
is on the order of .01. The probability that the total output of the 
transponder will exceed 16 dbm is of this same order (.01). Thus, the 
transponder is operating well within its designed values under congested 
conditions, and appears to be somewhat conservatively designed. 

In the receiver, a dynamic range of 45 db should be adequate to prevent 
saturation, with high probability. There is a low probability (4 x 10~ 5 ) 
that the noise plus the interference power will exceed a threshold level 
providing a 60 second warning at a 240 knot closing velocity. Thus, alarms 
caused by interfering signals are expected to be negligible under the 
conditions represented during the simulation. The cross«-product power 
term is the most significant interference source. 

Using the threshold level mentioned above, the system was in an alarm 
status 2.4 percent of the time. The probability that the second largest 

-4 

received signal will exceed the threshold level is on the order of 4 x 10 
under hour 11 conditions. While the receiver is an alarm status, the 
probability that the largest signal to noise plus interference ratio will 
exceed 20 db is on the order of .9. No cases were observed in which the 
signal to noise plus interference ratio was less than 10 db for the signal 
causing the alarm. 

Detailed study of the distribution of power levels from the aircraft 
in the population indicated that the power levels are approximately 
distributed normally. Actual distributions were found for the five closest 
aircraft to a randomly selected aircraft, and knowledge of the theoretical 
distributions permits extrapolation to higher density terminal area models. 


83 



In the flight tests, comparison of simulated and experimental values 
of range measurements indicated differences on the order of 10-15 percent 
over the calibration range of the equipment. This degree of accuracy 
permitted the determination of the sources of the measurement errors in 
the flight systems. Analyses of these sources of error were conducted 
by LRC personnel, using the results of the flight tests simulation. 

A comparison of idealized systems on the basis of warning times 
achieved for a given probability of alarm indicates that achieving the 
required protection without excessive alarms is an extremely difficult task. 
Systems based on proximity alone and on range and range rate measurements 
only are felt to be not feasible. As a minimum, altitude difference data 
should be used in the warning logic. 

Use of approximate miss distance discrimination is effective in 
reducing the number of false alarms, however, it was found that dangerous 
suppression of the warning time can occur under accelerating conditions, 
and that discrimination based on relative range acceleration should not be 
used in warning logic. The use of approximate miss distance measurements 
as assistance to the pilot in evasive maneuvers should be investigated 
further in future work. 

The extended modified tau criteria developed in Appendix A should have 

a beneficial effect on the false alarm performance, while not adversely 

effecting the system protection. This extended criterion, when simulated 

using the hour 11 data base, reduced the percent of time in alarm status 

from 2.2% to 1.4% for the criterion r <35 secs, and alt. diff. < 500 ft. 

m 

In summary, we can conclude that saturation and interference is not 
a problem with the system as presently designed, and under the flight 
conditions represented by the most congested hour of the data base. 
Simulations of the flight tests indicate that, in most cases, the sources of 
measurement error in the system are well known and can be accounted for 
satisfactorily. Comparisons of idealized systems indicate that the alarm 
threshold should be set to provide an alarm on the order of 30 seconds and 
that the modified tau (or extended modified tau) criterion should be used 
in warning logic, along with altitude discrimination on the order of ± 500 
feet. 
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APPENDIX A 

EXTENSIONS OF THE "MODIFIED TAU" WARNING CRITERIA 


Th'e "modified tau" criteria provides a warning based on the measure- 
ment of range (R) and range rate (R) , and an alarm is sounded if 

_ -R 4- [R 2 + 2 UR] ^J 2 

T m U < T mk (A— 1) 

where U and t ; are selected constants. The criteria, as shown in 
mk 

Appendix A, indicates that it may be. possible for the two aircraft to 
collide within x seconds if they make the worse possible maneuvers 
with a constant relative acceleration of U. 

This is a "worse case" criteria that requires relatively large 

9 

regions in the R, R plane to assure reasonable warning times x 

To obtain a more realistic idea of the warning time to be expected 
using this criteria, and assuming that pilot intent is unknown, it seems that 
probabilistic description should be used for the acceleration constraint 
U. That is, we consider the maximum relative acceleration as a random 
variable, and determine the expected or average value of the time-to- 
collision, again assuming that the aircrafts make the wprse possible 
maneuvers . 
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The value of this approach is that a more realistic estimate of the 




degree of hazard is obtained, and hence smaller regions in the R, R 
plane are required to provide sufficient protection* The fact that the 
worse possible maneuver is still assumed provides an additional margin 
of safety. 

We find the expected or average value of x^ by use of the 
expression . 


E {x } 
' m 


r co 

x (u) f(u) du 

J m 

— oo 


(A-2) 


where f(u) is the probability density function of maximum acceleration, 
considered now as a rahdom variable. t^Cu) is given by Equation A-l with 
U = u. As an example of the calculation, assume that any value of u is 

l 

equally likely, 0 < u < U, so that the probability density of u 
is 


f(u)=~, 0 < u < U 

The expected value of t is then 

m 



(A-3) 


(A-4) • 
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This integral is evaluated to give 


where 


T 

em 


E{x } = 
m 


2R 

D 


(4> (U) - 1) + log 


2 

<j>(U) + 1 


<KU> 




(A-5) 


and is used to designate the expected value of time-to-collision 

under the assumptions made. Curves of constant t have been plotted 

em 

in the R, R plane in Figure A-l for comparison with values of t . As 

m 

may be seen, the curves have the same general shape of the t curves 

m 

except for the region where closing velocity is small. For a given 


measured value of R and R, we find, as would be expected, that values of 


t are considerably less than value of x . 
em J m 

.The reliability of the estimated value x em could be improved if 
the assumption as to the density function of U could be removed. 


Further Extensions 

The modified tau criteria contains an unrealistic feature in that 
no distinction is made between normal and axial accelerations of the 
individual aircrafts involved. Axial (along path) acceleration components 
are known to be small in comparison to the normal acceleration components. 
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particularly in the terminal area with imposed speed limits. Also, one 
feels intuitively that hazards due to aircraft accelerations become less 
important as the measured closing velocity increases. That is, if a 
large closing velocity is measured, a near head-on encounter is 
indirectly implied. 

Thus, we are led to formulate mathematically a criteria that takes 
into account additional constraints on the allowable relative accelera- 
tion profile* 

Consider the diagram of Figure A-2. Here the present positions of 
the protected and intruding aircrafts are indicated by 0 and P 2 (T) 
respectively (i.e. aircraft 1 is at the origin of coordinates defined 
by unit vectors 1^ and 1^, with 1^ in the direction of p^CT)). The 
positions of the aircrafts at some time t later are given by 


T+t 

s 

T+t 

: 

m • 

' f 


P^Ct) dxds = 


J j 

. 

4 

T 

T 

T 


A^(x) dxds 


(A-to) 


and 


T+t s 

* t 

P 2 (t) dxds 

T T 


T+t' s 

# f 

A^Ct) dxds 

T T 


(A-7) 


where A (t ) and A_ (x ) are the acceleration functions of aircrafts 1 and 

1 * <4 

2, respectively. 
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Evaluations of the integrals on the left hand side gives 


P X (T + t) = 


T+t s 

„ /* f 

P X (T) t + 

J 4 

T T 


A^(t) dtds 


(A- 8) 


p 2 (t + t) 


P 2 (T) + p 2 (T) t + 


T+t s 

1 r 

A 9 (t) dxds 
J z 
T T 


(A-9) 


Now, since if a point can be reached in a given time by any sort of 
acceleration profile, it can also be reached by a constant acceleration 
of no greater magnitude, we assume a constant accleration to obtain 


P;l (T + t) - p x (T) t + A x t 2 /2 (A-10) 

P 2 (T + t) = P 2 (T) + P 2 CT) t + A 2 t 2 /2 (A-ll) 

These equations define the sets of points reachable by the aircrafts in 

a time t where A varies over all allowable constant vectors and t varies 
e 

from 0 to t . 

e 

For the set of allowable A^ and A 2 vectors, we assume that axial 
components are negligible and that the relative magnitude is limited to 
U, Thus the acceleration constraints are 

\ . 7 1 (t) = o 

A 2 • P 2 (T) = 0 (A-12) 

and |A 2 " A 1 | < U 
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For a collision to take place within a time t £ , we must have 

P 2 (T + t) = p^T + t) 


(A- 13) 


for some t, 0 < t < t . Using Equations A-10 and A-ll, we obtain 


r(T) + r (T) t + (A 2 - A 1 ) t 2 /2 = 0 

_ ± ± (A-14) 

where r(T) ■ p 2 (T) and r(T) = p 2 (T) - P 1 (T). 

Now resolve Equation A-14 into components normal to and along the 
relative range vector. 


l r Jr(T) + R(T)t - A 2 t 2 /2 sin & 2 - A^/2 sin 0^ = 0 
1 0 [*V n (T) t - A 2 t 2 /2 cos 0 2 + A^t 2 /2 cos 0^ 1 = 0 


(A-15) 


where, dropping the T notation for simplicity, R = 1 . r, R ** 1 . r, 

- - - - T r . A^ r 

V n = I 0 . r, A 2 = |a 2 |, A 1 = |A 1 |, sin 0 1 = and sin & 2 = 

T r ’ X 2 


Each vector component must equal zero, and squaring and adding the 


components gives 


(A-16) 


(R + Rt) 2 + (V t) 2 = 
n 


(A-, *+* ) A- 

- 4 - — ' " “V COS (0 2 + 9 1 ) 


t 4 . 
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The term in the brackets is, of course, the square of the 
magnitude of the relative acceleration vector. We desire to place 
additional constraints on the magnitude of this vector based on the 
fact that only components of acceleration normal to the velocity vectors 
are being considered. This appears impossible in general unless the 
relative headings of each aircraft is known, or unless additional 
assumptions are made. Following the latter course, we assume that 


P X (T) 


P 2 (T) 


(A-17) 


and that V^(T) =* 0. The first assumption is reasonable in the terminal 
area with imposed speed limits, while the latter assumption is 
equivalent to that assumed in the derivation of the "modified tau" 
criteria* 

With these assumptions, we note that 0^ = 0^ and that the maximum 
relative acceleration vector that will provide a solution is obtained 

when = A^. In reference [4], it is shown that if a solution to 
Equation A-16 exists for some ] A^ - A^| and some t, 0 < t < t g , it also 
exists for all greater values of l A 2 “ A^j. We have 

2 

/-I O A \ 

(A-18) 


A, 


max. effective |a^ ~ A.J < — (1 - cos 20) 


and since 2A^ < U if A^ = and |A^ - A.J < U, we have also 


S - IT 

(1 - cos 20) < (1 - cos 20) 


(A-19) 




* ' 
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Thus s under our assumption* a hazard will exist if 


(R + Rt) 2 



(1 - cos 20) 


(A-20) 


for some t, 0 < t < t . Taking the square root of both sides gives 


R + Rt = — — 
2 / 2 “ 


1/2 

[1 - cos 20] t 


(A-21) 


or 


R + Rt 



sin 0 


(A-22) 


The remaining task in developing this criteria is to estimate 
sin 6 from potentially measurable data. If we assume a reasonable 
maximum relative velocity exists, say V , (e.g. about 400 KTS in 
the terminal area) , then an estimate of sin 0 based on a range rate 
measurement is 


R ^ _ 

sin 0 s [1 - ( ^ ) ] R < V (A-23) 

max 


= 0 


R > V 


max 


Since cos 6 = — - (when |V.J = | V 2 1 and 0^ = 0^ as assumed), 
max ^ 

Thus, we have developed a criteria that has the desired behavior 
as initially discussed. At low values of R, the criteria behaves as 
"modified tau," while at large value of R, the criteria reduces to the 
"tau" criteria (t = R/R) , The criteria is summarized as: a hazard 

exists if there exists a t, 0 < t < t such that. 
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. (A-24) 


This criteria is plotted in Figure A-3 for several values of t 

e 

for comparison with the modified tau and tau criteria. 
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APPENDIX B 


THE D.C. OUTPUT OF A PRODUCT DETECTOR-LIMITER CIRCUIT 
WITH MULTIPLE INPUT SIGNALS 

Introduction 

The detector circuit used in the collision warning receiver has a 
configuration as shown in Figure B— 1. The Doppler tone signal is differenti- 
ated and split into two parts, one of which passes through a limiter and 
thence is combined with the other in the product detector. A low pass 
filter is used to eliminate high-order harmonics and smooth the output 
voltage. 



Figure B-l. Detection block diagram. 


Analysis 

Assume that four doppler signals of magnitudes P,Q,R, & S and noise 
n(t),are present at the doppler output, i.e. 


e(t) = P cos pt + Q cos qt + R cos rt + S cos st + N(t) 


(B-l) 
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where p,q,r, & s are the doppler frequencies. The output of the differ- 
entiator e^(t) is 

e^t) = P p sin pt + Q q sin qt + Rr sin rt + Ss sin st + n(t) 

4 

p,q,r or s <_ 2 tt x 10 


note that the derivative of the noise term is indicated only. Additional 
comments will be made in the subsequent analysis, for the case in which 
n(t) corresponds to white noise. 

The signal e^(t) is also the input to the limiter. The output of 
the limiter will now be derived, based on an extension of the results 
given by Rice [ 2 ]. The auto-correlation of the output from the limiter 


$ (T) is 


* l (t) 


CO CO 

• /* 

‘ J j f(e a > f(e b } P^a’V de a de b 

— OO —00 


(B-3) 


in which the output voltage of the nonlinear element is e ^ - f(e^) and 
p(e a e^) is the joint probability-density function of e2» that is 
e^(t) = e a (t + t) for the ergodic case. This equation will now be 
expressed in terms of the transform of f(e), which is defined as 


F(ju) = 


f(e)e ^ ue de 


— OO 


(B-4) 


Substituting (4) into (3) and interchanging the order of integration 
gives 


V T> ‘ ~2 j 


F(ju) du 


F(jv) dv 


J J 

—00 — OO 




jue a + jve b 


de de, 
a b 


(B-5) 
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(the constant term arises from inversion of (4). The last term in this 
equation may be recognized as the characteristic function of and v^, i.e., 


M(ju,jv,x) 


J 


P(e a e b ) 6 


jue a + jve b 


de de, „ 
a b 


(B-6) 


Thus (5) can be written as 
1 


» L ( T) - 2 

4t r 


F(ju) du 


F(jv) M(ju,jv,x) dv . 


(B-7) 


Assuming that the individual terms in (2) are independent, the charac- 
teristic function can be written in terms of the component sinusoidal 
and noise terms, 


M(ju,jv,x)= M p (ju,jv,x) MQ(ju,jv,T) I^Cju.jv,!) M g (ju,jv,x) M n (ju,jv,x). (B-8) 
Consider a typical term such as M r (ju,jv,x), by use of the ergodic hypothesis, 


M R (ju,jv,x) = T ^ i 


rT- 


exp[ju R cos rt + jv R cos r(t+x)] dt 


and by use of the relationship, 


u cos r t + v cos r (t+x) = /u 2 +v 2 +2uv cos rx fcos rx + 0) 


(where 6 is the phase angle), then integration yields 


Vju.jv.x) = J Q (Rr R ) 


(B-9) 


The notation F used in the above equation is defined by the expression, 


r = 2 2 

x u + v + 2uv cos xx 
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One other result is needed at this point. The characteristics func- 
tion for the gaussian noise voltage is given by. 


M 

u 


(ju. 


jv,x) - exp [ 


jfco 

2 


(u 2 + V 2 ) - ip uv] . 

T 


(B-10) 


In this equation is the correlation function of n(t) and i|> is the 
mean-square noise. The autocorrelation function, in terms of the above 
quantities, is 


P L <t) = ~^2 

L 4tt 2 


F(ju) e 


\bo 2 

* u *,f 


4>o 2 

- -2— v 


F(jv) e 


(B-ll) 


-i(i uv 

* 6 1 J 0 (P V W W J 0 (Sr B ) dV 


By use of the expansion 


Mr (ju, jv,t) = £ e (-l) n J (Ru) J (Rv) cos nrt, 

•tv- n n n 

n=o 

the autocorrelation function (11) can be written as 


00 oo oo 00 oo 


♦ L e> - I I M I kt £ n 

f — o g— o h=o i=o k-o K * 1 t 8 tllK 


cos f pT cos g qT cos hr T cos i sT. 


(B-12) 
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The term in this equation is given by 


. f +g+h+i+k f 00 

H = J 

fghik 2 it 


r 00 

j F(ju) u k J f (Pu) Jg(Qu) J^CRu) 


& u 2 

2 


(B-13) 


J ± (Su) e 


du. 


As Rice has pointed out, the dc and periodic terms in this expression are 
obtained by letting T 


OO 00 OO 00 


->• 00 " I l l l e n H fahi COS fpT c0s gqT 

f=o g=o h=o i=o 6 


x cos hr t cos i s t 


(B-14) 


This equation represents the autocorrelation function for the signal at 
the output of the limiter. Specific terms will be examined following 
analysis of the product detector. Note that the above expressions assumed 
an input of the same form as (1) 9 for an input such as (2) each of the 
P 9 Q, R, and S terms should be replaced by Qq s etc. 

The next task to be discussed is that of deriving an expression for 
the dc output of the product detector. The terms related to the dc output 
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can be identified by considering the following analysis * For a 
signal 


K- w, cos w,t 
Id d 

acting as an input to the product detector, and a square wave switching 
voltage 


CO 

OTr v sin tt/2 

2K 0 ) 77\ cos nw,t, 

2 L - n7T/2 d 

n=l 

the output will comprise the product of these terms. The only term con- 
tributing to the dc component is 

OTr Tr sin tt/ 2 r ^2 

2K 1 K 2 d tt/2 lcos nw d tl 

- 4 W 1 ' cos 2 ”d t)- 

The dc output of a product detector, in which the input and switching volt- 
age consists of incommensurable sinusoidal terms, is the Siam of fundamental 
coefficient products. The expressions for the harmonic and cross product 
(fluctuating) terms is extremely messy. These terms will not be considered 
since the present purpose is to evaluate distortion in the detected doppler 
signal arising from multiple targets. 
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The low-pass filter, shown in Figure B-l, uses the circuit shown below 


c 



The transfer function of this filter is 


e 7 (S) 


R 


_2 

sC 


R. 


2 sC 


R-. 



1 

sCR 2 + 1 


(B-15) 


The transfer function for a limiter may be found in Middleton p. (639) [9j, 


i 


.e. 


-uR 

1 - ° 

F(ju) = 23 — . 

2 

-u 

The symbols used in this equation are shown in the following sketch. 



Figure B-2. Definition of 3 
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Returning to (14) , the fundamental terms in the output of the limiter 
can be seen to result from equating one of the fghi indices to unity 
and the others to zero. Thus the magnitude, for the input sinusoidal 
P cos pt, of the output fundamental term is 


H 

P 


i 

7T 


2g 


1 - e 


-u 


-uR. 


° J L (pPu) J o (qQu) J Q (rRu) J q (sSu) 


du 


(B-16) 


2 

This equation also assumes that << R q . The contour for this integra- 
tion is shown below. 



The dc output of the product detector (using k e as the detector proportionality 
factor) is the product of the input signal and the coefficient of the 
fundamental terms, i.e,. 


E dc - k d + ^ H q + rRH r + sSH s ] ‘ (B-17) 

This equation represents the dc output of the sum channel. The outputs 

from the other channels are obtained by substituting the proper channel 

voltages for the P, Q, R, and S terms in the above equation (e.g,, page 9 
of ref [3]). 

Integration of Equation B-16 is accomplished as follows. 
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2H = — 

7T 


-iuR 

2(3 X - e — — J^) J q (Qu) J q (Ru) J q (Su) du 

-u 


F(PQRSU) 


_ 2ig 

7F 


r l-(cos u R - i sin uR ) 


F(P 9 Q, R 9 S, U) du 


”U 


u - a + j 6 


= ISi 

7T 




+ i 


r sin R o aF()da 
2 

-u 


= 0 since is odd 


- 0 since is odd 


r°° sin R u 

2H = J..(Pu)J (Qu) J (Ru) J (Su) du 

7T 2 1 O O O 

* u 


, q I sin R u 
4g o_ 

TT 2 

0 U 


J,(Pu) J (Qu) J (Ru) J (Su) du 
1 o o o 


(B-18) 


Note that the terms such as qQ have been abbreviated simply as Q in this 
integration. 

Numerical Calculations 

Computer programs have been written to solve the integral given in Equation 
B-18.- This integral calculates the fundamental output of the limiter, or 
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normalized gain for the particular input of interest. Curves have been 
calculated for various values of input signal level as shown in Figures B-3 
and B~4. Figure B-3 represents the case for two input signals and Figure B-4 
shows the normalized limiter gain for the case when a third signal is 
present. 

In the computer simulation, the dc voltage from the product detector 
is calculated using Equation B-17 , with Equation B-17 generalized to include 
all significant signals. 
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APPENDIX C 


PHASE PLANE PLOT OF A SET OF COLLISION TRAJECTORIES 
UNDER AN ACCELERATION CONSTRAINT 

It is informative, in investigation of the protection afforded by 
various systems, to visualize the relative trajectories in the phase 
(R, R) plane. 

For acceleration limited trajectories, we have 


|r| < U or -U < R < U (C-l) 

where R is range acceleration. Hence, the slopes of the trajectories in 
phase space are given by. 



U 

< a < *- 


(C-2) 


A set of allowable trajectories are sketched in Figure C-l, for the specific 
case of U = 1/2 g. (Corresponding to an aircraft bank angle of about 27° 
for one aircraft turning) 
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CLOSING VELOCITY* KTS 



0 1 2 3 4 5 

RELATIVE RANGE 9 NAUT. MI. 


Figure C-l. Phase plane plot indicating slope isoclines for trajectories 
with accelerations not exceeding ± 1/2 g. A set of 
allowable trajectories is also sketched. 
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APPENDIX D 

RECEIVER DATA PROCESSING 


General Description 

The following documents the calculations for voltages from the receiver 
that are analogs of range, range rate, range acceleration, normal velocity, 
tau, beta, miss distance, and modified tau. All of these voltage levels are 
calculated from the received power levels and frequencies. The calculations 
are based on both theoretical models and empirical data derived during the 
receiver calibrations at LRC. 

After calculation of the voltage analogs, a conversion is made to the 
corresponding geometrical value and the errors relative to the ground based 
radar values and experimental values are calculated. The inputs to this 
part of the simulation program include the received power levels and 
frequencies, the geometrical and experimental values of parameters, and the 
applicable system parameters. 

Input Parameters 

Inputs used in this section of the program include the following: 


At “ Calculation time increment 

- AGC loop time constant 

t,- - Discrimination time constant 

dl 


T 


d2 


= Differentiator time constant 


f(P,f^) = Receiver response characteristic as a 

fraction of IF amplifier input power 
level (P) and doppler frequency (f^) 
(experimentally determined) 


Experimental (determined from reduction of flight records) 
values of range, range rate, etc. 

Geometrical (derived from ground-based radar) values of 
range, range rate, etc. 
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Calculations 


Using the values of doppler frequency and power at the IF amplifier 
input* the voltage from the receiver-detector is calculated from the 
function f(P, f^), This experimentally determined receiver calibration 
function is shown in Fig, D-l. Figure D-2 indicates the detected voltage 
versus power at the IF amplifier input for a particular frequency (1600 Hz). 
The output of the detector is a voltage analogous to range in nautical 
miles. It should be noted that the function shown in Fig* D-l is 
introduced into the computer as points at 100 Hz intervals of frequency* 
and linear interpolation is used to calculate values at intermediate 
frequencies. 

The time response characteristics are introduced into the calculation 
by the following technique: The detector low-pass filter output voltage is 

given by the equation 



(D-l) 


where is the input voltage, x & is the filter time constant, and S is the 
Laplace variable. In differential equation form. 


T 

a 


dv 
o 

dt 


+ v 

o 


v i 


(D-2) 


Approximating the derivitive gives 


n n-1 
v - v 
o o 

T a At 




(D-3) 


where the superscrips indicate the calculation increment. Solving for the 
current (nth) value of output voltage gives 



(volts) 


CD-4) 
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Figure D-2. Power at IF amp input (DBM) * 
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Equation D-4 is then used to calculate the detector output (range) voltage. 
Provisions are, of course, made to assure the proper initial conditions to 
prevent transients. Similar expressions are used for the discriminator and 
differentiator filter voltage outputs. 

The percentage errors between measured, geometrical and experimental 
values (see printout of Fig. 18 ) are calculated as, for example, 

% error EXP-GEO - x 100 . (D-5) 

Geo. range 

Percentage errors due to power variations (e.g. gain calibration and 
antenna pattern characteristics) are calculated from expressions such as, 

% error (antenna pattern) = K (AG) (D-6) 

where AG is the antenna gain variation in db and K is a constant reflecting 
the detector slope (i.e. the "range" voltage per DBM power at the IF 
amplifier input) and a conversion from base 10 logs to base e logs. The 
derivation of Equation D-6 is as follows: The detector characteristics at 

1600 cps are given by 


log R = .02 P, + C, 
dbm 1 


(D-7) 


where R is the range in n. mi. and C is a constant. Differentiation gives 


dR 

R 


.02 

log 10 e 


dP 


dbm’ 


(D-8) 


and Equation D-6 follows immediately since a change in gain is equivalent 
to a change in power level. 
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The residual errors are calculated as indicated by the numbers on 
the printout sheet. For example, 


Residual (1 -2*3*4* 6) = 


' (EXP-GEO % error + 100) - 

100 


(dop. fil. % error + 100) 
100 


(det. slope % error + 100) 
100 


(el. pattern % error 4- 100) 
100 


x 


(gain calib. % error + 100) 1 ion-100 

100 J 
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APPENDIX E 
SYSTEM PARAMETERS 


System parameters form an input to the simulation program, and are 
a punched card input. Following is a list of the system parameters used. 


No. 

Parameter 

Units 

Value 

1 

Transmitted Power 

DBM 

30.00 

2 

High Transmitter Freq. 

MC 

4252.50 

3 

Low Transmitter Freq. 

MG 

2702.50 

4 

Oscillator Deviation Del 

MC 

2.50 

5 

Oscillator Deviation Del 

MC 

.50 

6 

Bal. Mixer Bias Power 

DBM 

-22.00 

7 

Transponder Output Sat. Power 

DBM 

30.00 

8 

Transponder Gain Constant . 

DB 

176.50 

9 

Transponder Input to Bal. Mixer Gain 

DB 

49.50 

10 

Transponder Input Channel Bandwidth 

MC 

12.00 

11 

Transponder Output Channel Bandwidth 

MC 

2.00 

12 

Transponder Noise Figure 

DB 

8.50 

13 

Transmitter-Transponder Isolation 

DB 

-150.00 

14 

Receiver Noise Figure - 

DB 

6.30 

15 

Doppler Filter Bandwidth 

KC 

3,00 

16 

Power Gain Constant-Rec . Input to I.F. 
Amp Input 

DB 

28.50 

17 

Transponder-Receiver Isolation 

DB 

-150.00 

18 

Spare 

— 

0.00 

19 

Spare 

— 

0.00 

20 

AGC Time Constant 

SECS 

3.00 

21 

Discriminator Time Constant 

SECS 

1.00 

22 

Differentiator Time Constant 

SECS 

1.00 

23 

Spare 

— 

0.00 
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No, 

Parameter 


Units 

Value 

24 

4,25 Transmitter Antenna 

Gain 

DB 

3.90 

25 

4.25 Transponder Antenna 

Gain 

DB 

-1.50 

26 

2.70 Transmitter Antenna 

Gain 

DB 

5.20 

27 

2.70 Transponder Antenna 

Gain 

DB 

-1.00 

28 

1*55 Receiver Antenna Gain 

DB 

1,80 

29 

1,55 Transponder Antenna 

Gain 

DB 

.50 

30 

Spare 


— 

0.00 


* 

Note i Antenna gain values are adjusted for loss factors in some cases. 
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APPENDIX F 
FLIGHT TEST DATA 


In this Appendix, computer printouts for the flight tests discussed 

in Section V are presented. Although the printouts were obtained every 

5 seconds during each flight test , the printouts are given only for the 

first point at which the alarm was received in order to conserve space. 

In cases where no alarm was received (in the simulated flight) a printout 

2 

near the point of minimum R t is given. For interpretation of the 
printouts, see the discussion in Section V-B. 

Is should be noted that during the flight test, the systems in the 
aircrafts were slightly different. For simulation purposes, however, 
both systems were assumed to be configured as was the prime system. Thus 
the simulated systems (system 1 and system 2) are identical. System 1 
was the prime test aircraft, and the aircraft on which the experimental 
data were obtained. 
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Figure F-l. Computer printout o£ flight test simulation for Flight 7A, 
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Figure F-2. Computer printout of flight test simulation for Flight 7B. 
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Figure F-3. Computer printout of flight test simulation for Flight 7C 
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Figure F-4. Computer printout of flight test simulation for Flight 7D 
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Figure F-6. Computer printout of flight test simulation for Flight 8C. 




Figure F-7. Computer printout of flight test simulation for Flight 9B. 
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Figure F-9. Computer printout of flight test simulation for Flight 9D. 
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Figure F-10. Computer printout of flight test simulation for Flight 9E. 
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Figure F-12. Computer printout of flight test simulation for Flight 10B. 
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Figure F-14. Computer printout of flight test simulation for Flight 1GD 
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APPENDIX G 

ANTENNA PATTERN MODELS 
Pattern Models for Statistical Studies 

For the statistical studies, the actual antenna patterns supplied 
RTI by NASA-LRC were represented by analytical expressions* Figures 1 
through 7 show the static patterns with' a broken line superimposed to 
indicate the corresponding analytical expression* Equations 1 through 9 
below are the corresponding analytical expressions (note that the 
expressions indicate absolute voltage gain and are converted to power 
gain in db in the program) * 


Transmit ter -Azimuth 
(4252.5 & 2702.5 MHz) 

. _ 1 -0.479(130 - 6) 0,314 ; 0 < 6 < 130° 

“ e 

A = 0.039 ; 130° < 0 < 180° 


Tr ansmi t ter-Az imuth 
(4252.5 & 2702.5 MHz) 


sin 96 
0.162 0 


; 0 <_ 0 <_ 1° 

; 1° < 0 < 180° 


Transponder-Elevation 
(4252.5 MHz) 


A = cos 3.83 0 


; 0 < 0 < 10° 


sin 4.87 8 
0.00972 0 


; 10° < 0 < 90° 


(G-l) 


(G-2) 


(G-3) 
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Transponder-Elevation 
(2702.5 MHz) 


A = cos 3.74 0 

. sin 40 

^ 1 43 

0.0301 6 J 


Transponder-Elevation 
(1550 MHz) 


A = 1 

A = sin 46 
0.0722 0 


; 0 £ 6 £ 10 ° 

; 10° < 0 < 90° 


; 0 <_ 6 <_ 1 ° 

; 1° < 6 < 90° 


Receiver-Azimuth 

A _ x _ e -O.O371(156-0) 0 * 847 

A = 0.056 


; 0 < S < 156 

; 156° < 0 < 180° 


Receiver-Elevation 

A = cos 2.34 0 

. sin 3.756 

A TTa 

0.00773 0 ’ 0 


; 0 <_ 6 <_ 20 ° 

; 20° < 0 < 180° 


(G-4) 


(G-5) 


(G-6) 


(G-7) 

(G-8) 


Pattern Models for Flight Test Simulation 

This section describes the antenna pattern models used for the flight 
test analysis. The horizontal and vertical plane patterns were subjected 
to a Fourier analysis to determine the Fourier coefficients for each of the 
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patterns supplied RTI by NASA-LRC, The gain at a desired aspect (0,<j>) 
was then reconstructed from the series representation (ref. G-l), 


G( e ^) db " G( 6 ) db + G ( ) db 


20 



p=0 


20 

+ Z 


P~0 



COS 


2^6p 0 

2N+1 + % 



cos 


, A 
2N+I + 3 p 


sln iHi] 

3ln l»i] 


(G-9) 


where , 

0 is azimuth in degrees (0 <_ 9 <_ 360) pos. CW 

<j> is elevation in degrees (0 <_ 8 360) pos. up 

from nose 

a ,6 are the fourier coefficients 
P P 

N is no. of data points/2 (=180) 

A brief analysis was conducted to indicate the total number of harmonics 
required to accurately describe the pattern signature. It was observed that 
10 harmonics did not adequately describe the null regions while it was felt 
that in excess of 20 harmonic would introduce computing errors of sufficient 
extent to actually degrade the pattern. For the purpose of flight test 
analysis, 20 harmonics were used in the calculation of the gains. 

The reconstructed patterns (see Fig. G-8 for an example) were observed 
to reproduce the static patterns within approximately two db over the major 
portion of the lobing structure and suffered some degradation in the null 
regions due to lack of high frequency components. This accuracy is considered 
adequate for the analysis 


Ref. G-l Goertzel, G. , "Fourier Analysis," Mathematical Methods for 
Digital Computers , Wiley 1964 
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Az . 4252.5 MHz 
Transmit 



Figure G-l. Transmitter azimuth pattern (4*252 GHz) • 
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21 


Figure G-3* Transponder receiver * elevation (4.252 GHz) 







El. 2702.5 MHz 

Transponder 

Receive 



Figure G-4. Transponder receiver 3 elevation (2*702 GHz) 
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Figure G-5* Transponder transmitter;, elevation (1*5 GHz) 
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Figure G-6. Receiver elevation (1*5 GHz), 
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Figure G-8. Example of Fourier fit* 




APPENDIX H 


ANTENNA GAIN VARIATION WITH RELATIVE 
ELEVATION ANGLE 

To determine the variation of measured range with elevation angle, 
the antenna patterns obtained in the LRC anechoic chamber have been used 
to derive the curve of Figure H-l. This curve shows the ratio of measured 
range to actual range as a function of the elevation angle of the 
intruder. The curve takes into account the gain variation with elevation 
of all antennas in the system. 

As may be seen from the curve, the ratio of measured range to actual 
range remains relatively constant over a ± 10 degree interval, and then 
rapidly deteriorates. 
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RANGE MEASURED 
RANGE ACTUAL 





